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With the recent measurement of 2/^(10610) and 2^(10650) — > BB* +c.c. and B*B*, we investigate the transi- 
tions from the 2^(10610) and Z/,( 10650) to bottomonium states with emission of a pion via intermediate B B* 
meson loops. The experimental data can be reproduced in this approach with a commonly accepted range of 
values for the form factor cutoff parameter a. The T(35)7t decay channels appear to experience obvious thresh- 
old effects which can be understood by the property of the loop integrals. By investigating the a-dependence 
of partial decay widths and ratios between different decay channels, we show that the intermediate B B* meson 
loops are crucial for driving the transitions of Zi,/Z'i^ Y{nS)TZ with n= 1,2,3, and h/, {mP)lt with m = \ and 2. 

PACS numbers: 13.25.Gv,14.40.Pq,13.75.Lb 



I. INTRODUCTION 

Recently, two charged bottomonium-like structures Z,^(10610) and Z^*( 10650) (abbreviated to Z^ and Z'f^ in 
the following) were observed by the Belle Collaboration in the 7t*T(«5) (n = 1,2,3) and n'^hb{mP) (m =1,2) 
invariant mass spectra of T(55') T{nS)n'^n^ and ht,{tiiP)%^n^ decays Ullzll. The reported masses and widths 
of the two resonances are M7+ = 10607.2 ± 2.0 Me V, Ty+ ^ 18.4 ±2.4 MeV andM7,+ = 10652.2 ±1.5 MeV, 

■ ^b 



— 11.5±2.2MeV \mM- Analyses of the charged pion angular distributions favor the quantum numbers of the 
Z-states I^{J^) = 1 + (1 + ). Evidence for the charge neutral partner Z^ is found in a Dalitz plot analysis of T(55) 
T{2S)n\^ with 4.9a significance by Belle Collaboration Its measured mass M^o = 10609+g±6 MeV is also 
consistent with that measured in the charged mode. Since Z^'s are isotriplet states, they need at least four quarks as 
minimal constituents, which makes them ideal candidates for exotic hadrons beyond the conventional qq mesons. 
Note that the decay rates of T(55') — > Z^K — > Y{nS)KK are comparable to those of T(55') Zi,K — s- lii,{mP)KK. This 
implies unusual dynamic mechanisms undergoing the decay process since the transition to lii,{mP) would require 
the flip of heavy quark spin and should be suppressed in the heavy quark mass limit. 

Before the observation of Z^ and Z^+, the authors predicted the existence of loosely bound 5- wave BB* molec- 
ular states lUIlt]. In Ref. the authors predicts the possible existence of Z?'*'/?^*' molecular candidates within 
one-boson-exchange model. Since the Z^ and Z,'^ are charged and close to the BB* and B*B* thresholds, many 



studies show that they could be 5-wave BB* and B*B* molecular states ll8l- ll4ll . In Ref. lll5ll . the masses of 5-wave 
heavy tetraquarks bubd and bdbii with 7^ = 1+ are extracted by the chromomagnetic interaction Hamiltonian, 
which turn out to be compatible with the corresponding masses of Z^ and Z'l^ . The QCD sum rule calculations 
provide a tetraquark interpretation f[l6ll . Meanwhile, the tetraquark picture is applied to the understanding of the 
decays of Z±/Zf 7i±T(«5) and K^hh{mP) |fl7|. 

Besides the spectrum study, the production and decay of Z^ and Z'l^ are also investigated extensively. Con- 
sidering Z'^ and Z'l^ to be BB* and B*B* molecular states, Voloshin estimates their production in the radiative 



decay of Y(5 5) fllSll . and t he p ion-emission transitions from Z'^ and Z'l^ to lighter bottomonia are investigated 
effective field theory assuming that Zf and Z(^ are the BB* and B*B* molecular states. 



by Refs. Ill8lll9ll . In Ref. Ill4tl . the properties of Zy and Z^ were studied in the framework of a nonrelativistic 



'Electronic address: | gli@mail.qfnu. edu.cn | 
^Electronic address: zhaoq@ihep.ac.cn] 



2 



The intermediate meson loop transitions have been one of the important nonperturbative transition mechanisms 
in many processes, and their impact on the heavy quarkonium transitions, sometimes called coupled-channel ef- 
fects, has been noticed for a long time |23-22]]- By applying the on-shell approximation, the bottom meson loops 
were s ugg ested to p lay an important role in the Y(55') transitions to the lower T states with the emission of two 
pions 12311 or one r\ 12411 . This mechanism seems to explain many unusual properties that make the Y(55) different 
from T(45). Similar approach was also applied to the study of Z/, and ZjJ, by Liu et al. ill ill . 

In this work, we will investigate the decays of Zb T{nS)K and Zh hb{mP)n via intermediate S-meson loops 
in an effective Lagrangian approach (ELA) with the favored quantum numbers I^{J^'^) = 1 + (1+^) for the Z^IZ'^^. 
We try to enhance the scenario by quantitative calculations that the bottomed meson loops are crucial for explaining 
the experimental results for and — > BB* + c.c. and B*B* , and Zj, and Z^ Y(n5)7t and h[,{mP)n. 

The paper is organized as follows. In Sec. |II] we will introduce the formulae for the ELA. In Sec. |III] the 
numerical results are presented. The Summary will be given in Sec.lIVI 





FIG. 1: The hadron-level diagrams for — > Tlt+ with s(*)b(*) as the intermediate states. Similar diagrams for and Zj^ 
states decays. 




FIG. 2: The hadron-level diagrams for — > /ii,7t+ with s(*)s(*) as the intermediate states. Similar diagrams for Z^ and Zj] 
states decays. 



II. TRANSITION AMPLITUDE 

In order to calculate the leading contributions from the bottomed meson loops, we need the lea ding order effec- 



tive Lagrangians for the couplings. Based on the heavy quark symmetry and chiral symmetry Il25ll26ll . the relevant 
effective Lagrangians used in this work are as follows, 

^inS)B(*)BM = i8TBBT^,{d^'BB-Bd^'B) -grB*B£"'"'%Ty{daB*pB + BdaB*p) 

-igTB'B* {r>'{d,B*^B* - B*^d,B*) + {d,T,B*^ - T,d,B*^)B*'^ + 3^^; - d,rB;)}ii) 

\(,nP)B(')BM = gh,B'Bhl{BBl + B^B) + igh,B'B'^'"'%J^bvBlBl , (2) 
^B'Bi'h - ~igB'Bn{B,d,P,jBj - B*%P,jBj) + \gB*B*n£'""'^B*^d,P,j'^aB%, (3) 

where = ) andB^*)^ = (b^*^',B^*^'^] correspond to the bottom meson isodoublets. 
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With the ex£erimental data for BR{Z+ B+B*° + B^B*+) = (86.0±3.6)% andBR{Z'+ ^ B*+B*^) = (73.4± 
7.0)% from f27'], we obtain gz,,B*B = 13.39 GeV and gz'B*B* ~ 0.32. The relations 



gZtB'B = -gZ,,B*B* mZf, 1 / , gz!. B'B^ Sz'. B'B* '"zl \ / , (5) 

are applied to extract the couplings for gz,jB*B* and gz'^B'B- 

In Eq. ([T]), The following couplings are adopted in the numerical calculations, 



I STBB ms* 

gTBB = lg2\/m^mB , gTB'B = , , gTB'B* = gTB'B \ '"B* , (o) 

y/mBinB* V '"B 

where g2 = \/ntx / {IniBfj)', fy and my denote the decay constant and mass of T(«5'), respectively. The decay 
constant /x can be extracted in T(«5') — > e+e^: 

nrinS)^e-e-) = '-^^, (7) 
27 OTr(«s) 

where a^u = 1/137 is the fine-structure constant. By adopting the mass values in Tableland data for the leptonic 
decay widths of T(«5) states: r(T(15) ^ e+e") = 1.340±0.018 keV, T{T{2S) e+e") = 0.612±0.011 keV, 
r(Y(35) ^ e+e-) =0.443 ±0.008 keV we obtain fy^,,) =715.2 MeV, /t(25) =497.5 MeV, and /t(35) = 
430.2 MeV. 

In addition, the coupling constants in Eq. Q are determined as 

niB* 

gh,,BB* = -Igx^mh^mBmB" , ghi,B*B*=2gi , (8) 



with gi — — V'^fco/^/Aio' where m^^^g and /j,^^, are the mass and decay constant of Xboi^P), respectively 1291], i.e. 
fx„o = 175±55MeV M, f^^i2P)/fx,„iiP) = /r(25)//r(i.), and/,^„(2^) = 121.6 MeV. 
The coupling constants relevant to the pion interactions in Eq. are 

2g , gB*Bn 

gB'Bn = —^JmBmB" , gB*B*K = , , (9) 

fit y/mBniB* 



where g = 0.44± 0.03;o:oo MM and /n = 132 MeV are adopted in this work. 

The loop transition amplitudes for the transitions in Figs. [T]and|2]can be expressed in a general form in the 
effective Lagrangian approach as follows: 

where (; = 1,2,3) are the vertex functions; a,- = qj — mf (; = 1,2,3) are the denominators of the intermediate 
meson propagators. We adopt the form factor, Yli !J'i{nii,qj), which is a product of monopole form factors for each 
internal mesons, i.e. 

X\ji{mi,q}) = 7i{mx,q\)72{m2.q2)73{mi,ql) , (11) 
i 

with 

where A, = m, + aAgco and the QCD energy scale Aqcd = 220 MeV. This form factor is supposed to parameterize 
the non-local effects of the vertex functions and remove the loop integral divergence. 



The explicit transition amplitudes for (p,) -> b'*' {qi)B^*\qi (^2)] — > T{nS) {p f)n{pTi) via those triangle 
loops are given as follows: 



I I 



2 2 ' n^K'g") (13) 



2 2 2 2 

q'Y — wjj — — m 



^BB'[B'] = (0 j ^^^fez6B*Be,>.]feT(n5)B*Bepa5T/'/e/°g'2][-^B*B*7tee,^K;tP^g'2] 



2 2 2 2 2 7 

Mb*B[B'\ = (O'^y^ ^^bz6B*Be/;,]feT(n5)B*B*(^pag4x-gpT^a4 + ^p4^ax)e7(?l +^2)^][-gB*B7cP7ce] 

,- ^^^^^^ ,,,, 

^j- — m| ~ ^"2 ^3 ~ '"3 I 

A^B*B' [B] = (0^ j \SZ,,B-B' e^vap^ [gl (nS)B' B^pa'qzP^f^Y q{\ [gB'B%P%(i\ 

i --(-^7A!/"^)na(».,..?) (16, 



q\ — m\ q2 — nij ^3 — '«3 



A/b*B*[B*] = {if J ^^^kz,,B'B'e^vapPfe)']fer(«S)B*B*fepa§5x-§px§o5+gp4gox)e7(?l +?2)^][-§B'B*7cee(t)KXPn92] 



/(-g"^ + g?gf/m^) ij-g'^' + q^.q'./ml) ij-g^^ + q^qt/mj) 



2 2 

q2 — q^ 



Also, the explicit transition amplitudes for Zi,{pi) B^*\qi)B''*\q3)[B^*\q2)] hh{mP){pf)K{pn) via those 
triangle loops are given as follows: 



^BB*[B*1 = {i? I 7KZ^lSZbB*B£i^,][ghk{mP)B*B£*f0][-gB*B*neQ^KXPKq2] 

-2 2 2 2 2 2 n^('"''^^) (1^) 



[gZi,B*Be,>,] [ghh{mP)B*B^*fp\ [-gB*B*n£e,^KkPn<l2] 

i K-g^' + qU/ml) K-g^'^ + ctAl^l) n rr^ 2^ 



q^ ^2 ~ ^2 ?3 



^B'B[B'] - {i)^ j Jj^[gZ,,B*BS-i^][ght{mP)B*B*^polzP^f^Y][-gB*BnP%%] 

q\ — m J g'j ~ ^"2 ^3 ~ '"s i 

(i'^^2 



Mb''B*[B] = (0 y ^^^fez,,B*B*e^va|5/5,eo]fe/!4(m/')B*Be/p]feB'Bit/5ji9] 



X 



i(-.-p+tf,?/„.;) ^ . .(-,^Y,g|/„.j) (20) 

<7j — OTj ^5 — ^3 — OT3 ; 

Mb*b'[B'] = (i)^ J J^^[gZi,B*B*£tjvapp'i£l][gh^{mP)B*B*^pai,TP^f^f"][-gB 
ii-g"^ + q'fqV>n\) i{-g-^ + q\q\lm\) /(-gP* + q\q\lm\) 



2 2 2 2 2 7 

O'l — m\ q^ — q\ — 'wj 



W7,(mi,q]), (21) 
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FIG. 3: (a) The a-dependence of the branching ratios of Z,+ T( 15)71+ (solid line), T(25)7l+ (dashed line) and T(35)7l 
(dotted line), (b) The a-dependence of the branching ratios of Y(lS)ll+ (solid line), Y(2S)ll+ (dashed line) and T(35)7l 

(dotted line). The experimental values are taken from ||27|| as a reference. 




0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 

a a 

FIG. 4: (a) The a-dependence of the branching ratios of Z'^ — > hh{lP)Tt'^ (solid line) and /?i(2/')ir+ (dashed line), (b) The 
a-dependence of the branching ratios of /?i,(lP)7t+ (solid line) and /ji,(2P)ll+ (dashed line). The experimental values 

are taken from ||27|| as a reference. 



where pi, p /, p^ are the four-vector momenta of the initial Zt,, final state bottomonium and pion, respectively, and 
qi , q2, and qj, are the four-vector momenta of the intermediate bottomed mesons as defined in Figs. [T] and |2] 



III. RESULTS 



TABLE I: A summary of meson masses adopted in the calculation. 

States T(15) T(25) T(35) hb{lP) hb{2P) B B* n 
Mass (MeV) [2^] 9460 10023 10355 9898 10259 5279 5325 140 



Proceeding to the numerical results, we list the meson masses involved in the hidden-bottom decays of Z^/Z^^ 
in Table |l] Several points concerning the determination of the form factor cutoff parameter which would be the 
only free parameter in a decay channel, should be clarified. First, we determine the cutoff parameter a for each 
channel separately by the experimental data. As shown in Table [III it is possible to find an appropriate range 
of a values for each decay channels that can account for the data via the intermediate bottomed meson loops. 
Meanwhile, one notices that the a values for /Z'^ — > T(35')7i+ are much smaller than other channels which 
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TABLE II: List of branching fractions for the and X'^ decays. The last column values are obtained at the average of the 
central a values exclude the T(35)7r+ channels. The experimental values are taken from ||27[l as a reference. 



Initial states 


Final states 


Exp. % 


a range 


a value 


This work % 


Z+(10610) 


Y(lS)it+ 
Y(25)jr+ 
Y(35)jr+ 

hb{2P)Tt+ 


0.32±0.09 
4.38± 1.21 
2.15±0.56 
2.81± 1.10 
4.34±2.07 


1 47+0.18 
^•^'-0.20 

1 76+''-28 


"•3^-0.09 
^•'°-0.30 

9 Qn+0.60 

^•^"-0.70 


1.81 


0.50 
4.57 
9.53 
3.03 
1.36 


Z^+( 10650) 


Y(lS)lt+ 
Y(2S)jt+ 
Y(3S)jt+ 

h{2P)%+ 


0.24 ±0.07 
2.40±0.63 
1.64 ±0.40 
7.43 ±2.70 
14.8±6.22 


1 9^+0.16 
^•^-'-0.18 
1 90+0.18 
^•^=^-0.21 

19+°"' 
"•^^-0.01 

1 36+°-2° 
l.JO_0 24 

1 62+°-^'* 
i -0^-0.43 


1.38 


0.31 
2.71 
5.55 
7.72 
11.18 



TABLE III: The branching ratios of decay rates for Z+ Y(35)7l+ and Z'+ Y(35)7l+ with Mb- =Mb = 5279 MeV (Scheme- 
I) and Mb' =Mb= 5325 MeV (Scheme-II). 





channels 


Z+(10610) ^Y(35)n+ 


z;+ ( 10650) ^Y(35)7t+ 


Scheme-I 


a 


T QA+0.76 
•'•°°-0.96 


■X 19+0.52 
■^•^^-0.58 


BR (%) 




1 64+0-40 
^•"^-0.40 


Scheme-II 


a 


-=■•^^-0.60 


T 08+0-49 
•'•"^-0.50 


BR (%) 




1 64+°-"*" 



indicates some unusual feature with this channel. In Table HH we also list the a values for each decay channels that 
can reproduce the experimental data. An alternative test is that we make an average of the a values for the and 
decays separately without including the T(35)7t channel, and then check whether it is possible to describe the 
experimental data with single values of a for Zh and Z^, respectively. Interestingly, as shown by the sixth column 
of Table HIl with a = 1.81 and 1.38 for the Zi, and Z^ decays, respectively, the data can be reasonably accounted 
for except for the T(35')7i channel. 




1.0 1.5 




FIG. 5: (a) The a-dependence of the ratios R^^ (solid line), R|2 (dashed line), R^? (dotted line), R^? (dash-dotted line) defined 
in Eq. l l22t . (b) The a-dependence of the ratios rz^ (solid line), r^' (dashed line) defined in Eq. l l23t . 
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We also check the a-dependence of the decay branching ratios in order to give a quantitative estimate of the 
cutoff uncertainties in the loop integrals. The numerical results are summarized in Figs. |3]and|4]for the Zb and 
decays into T(/i5')7t and ht,{mP)K, respectively. 

In Fig.[3](a), we plot the a dependence of the branching ratios of T(15')7t+ (soUd Une), Y(25)7t+ (dashed 

line), and T(35')7i+ (dotted line), respectively. A predominant feature is that the a dependence of the branching 
ratios are quite stable, which indicates a reasonable cutoff of the ultraviolet (UV) contributions by the empirical 
form factor As shown in this figure, at the same a, the intermediate B-meson loop effects turn out to be more 
important in Z^ T(35')7i than in Z^ — )■ Y(15, 25)71+. As a result, a smaller value of a is favored in Z^ — ^ 
T(35')7t+. This is understandable since the mass of T(35') is closer to the thresholds of BB* or B*B* than the other 



two states 12811 . Thus, it gives rise to important threshold effects in Z^ — > T(35')7t+. 

One notices that the a-dependence of the branching ratios for Z^/Z^+ — > T(35')7t+ is stabler than those for 
T(15', 2S)n. This indicates that the enhanced branching ratios are not from the off-shell part of the loop integrals. 
As we know that the dispersive contributions become rather model-dependent near threshold, the enhanced (but 
rather stable in terms of a) branching ratios for Z'^ /Z'^^ — > T(35')7t+ suggests that more stringent dynamic con- 
straints are presumably needed to describe the near-threshold phenomena where the local quark-hadron duality has 
been apparently violated. What makes this process different from e.g. hcT\P in Ref. 13311 is that there is no can- 



celations between the charged and neutral meson loops. As a consequence, the subleading terms in Refs. 11331134 1 



become actually leading contributions. In this sense, a new power counting in the nonrelativistic effective field 



theory should be exploited for the Z/,/Z^ decays ll35ll 



Figure[3](b) presents the branching ratios of Z'j^ — > T(n5)Jt+, and the notation are the same as Fig.[3ja). 
The a dependence of the branching ratios of Z'^ jz!^ hi,{lP)K^ (solid line) and /!/,(2P)7i+ (dashed line) is 
presented in Fig. |4] The experimental data are denoted by points for corresponding decay channels. The data 
for Z^ hi,{lP)n+ and hh{2P)%+ can be reproduced with a = 1.76+o3o ^.QO+o^q, respectively. For Z^ 
hh{lP)z^ and hh{2P)%^, the values of a — 1.36^q24 ™d 1 -62^0 43 can be determined by the experimental data. 
As shown in Table Ull the decay channels for both Zb and Z^ hh{\P, 2P)k can be reasonably accounted for by 
the averaged values of a 1.81 and 1.38, respectively. Moreover, as shown in Fig.|4] the a dependence turns out 
to be stable for both Z/, and Z^ decays. The stabler behaviors for Z'^ and Z'/^ hh{2P)Ti^ than hh{\P)Ti^ indicates 
the closeness of the B*B* threshold to the /!/,(2P)7i+ threshold and the dominance of the meson loop contributions 
due to the open threshold effects. 

It would be interesting to further clarify the uncertainties arising from the introduction of form factors by study- 
ing the a dependence of the ratios between different partial decay widths. For the decays of Z^/Z^+ — > T(n5')Jt+, 
we define the following ratios to the partial decay widths of Z'^ jz!^ T(25')7t+: 

12 ^ r(Z+^T(15)7r+) 32 ^ r(Z+^Y(35)7t+) 
r(Z+^T(25)7t+)' r(Z+ ^ T(25)7t+) ' 

,2_ r(z;+^T(15)7r+) 32_ r(z;+^Y(35)7r+) 

r(z;+^T(25)7t+)' 4"r(z;+^T(25)7t+)' ^ ^ 



which are plotted in Fig.|5](a). The stabilities of the ratios in terms of a indicate a reasonably controlled cutoff for 
each channels by the form factor 

For the decays of Z^ jz!^ — > hbi^P, 2P)7i+, the following ratios are defined: 

^ r{Z+ ^ hb{2P)n+) ^ r{Z'+ ^ hb{2P)K+) 
TiZ+ ^ hbmn+y T{Z'+ ^ hbilP)n+y 

The a dependence is then plotted in Fig. |5](b), which also appears to be highly stable. Since the first coupling 
vertices are the same for those decay channels when taking the ratio, the stability of the ratios suggests that the 
transitions of ZbjZ'^^ Y{nS)K and hb{mP)% are largely driven by the open threshold effects via the intermediate 
B meson loops. In order to understand this, the following analysis is carried out. First, one notices that we have 
adopted the couplings for the hb and T to BB* otB*B* in the heavy quark approximation. Since the physical masses 
for B and B* are adopted in the loop integrals, the form factor will introduce unphysical pole contributions of 
which the interferences with the nearby physical poles would lead to model-dependent uncertainties. By assuming 
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Mb- =Mb= 5279 MeV (Scheme-I) and Mg* = Mb = 5325 MeV (Scheme-II), namely, by partially restoring the 
local quark-hadron duality, we calculate the partial widths of Z'^ and T(35')7t. We expect that the partial 

restoration of the local quark-hadron duality will significantly lower the partial widths since there will be only one 
physical pole in the loop and the unphysical one can be easily isolated away from the physical one. As a result, the 
inferences caused by the closeness of the unphysical pole will be reduced js^. As listed in Table HIH in the heavy 
quark limit, i.e. Mb* = Mb, the partial widths of and — > T(35')7t can be reproduced at much larger a. This 
is a rather direct demonstration of the sensitivity of the meson loop behaviors when close to open threshold and 
when the dispersive part becomes dominant. 

In brief, we find it is possible that with the same values of a for different decay channels, experimental data 
for the Z), and ZjJ, hadronic decays can be accounted for in terms of intermediate B meson loops except for the 
T(35')7t channel where the close-to-threshold effect plays an important role. Recognizing this is helpful for us to 
understand the experimental results, and establish the the intermediate B meson loops as the dominant transition 
mechanisms for the Zt, and Z^ decays. 

IV. SUMMARY 

In this work, we investigate hidden-bottom decays of the newly discovered resonances Z^ and Z'j^ via inter- 
mediate Z?-meson loops. In this calculation, the quantum numbers of the neutral partners of these two resonances 
are fixed to be l'^{J^^) = 1 + (1+^), which is currently favored by the experimental analysis. In the ELA, the 
present experimental data can be reproduced with a commonly accepted range of values for the cutoff parameter a 
except for the Y(35')7t channel where the close-to-threshold effect plays an important role in the process of Z^ and 
Z',,^B*B*{B)^T{3S)n. 

Our results show that the a dependence of the branching ratios are quite stable, which indicate the dominant 
mechanism driven by the intermediate meson loops with a fairly well control of the UV contributions. We also 
pointed out that the results become sensitive to the meson loop contributions when the final state mass threshold 
are close to the intermediate meson thresholds in our calculation. Namely, the effects from the unphysical pole 
introduced by the form factors would interfere with the nearby physical poles from the internal propagators and 
lead to model-dependent uncertainties. It is also a consequence of the violation of local quark-hadron duality. 
Such a phenomenon has been discussed in Ref. 1B4I1 . Further experimental and theoretical studies llssll of the 
consequences from such an intermediate meson loop effects would be important for providing more quantitative 
information on the structure of Z'^ and Z'/^ . 
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